Index of refraction of molecular nitrogen for sodium matter waves 
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We calculate the index of refraction of sodium matter waves propagating through a gas of nitrogen 
molecules. We use a recent ab initio potential for the ground state of the NaN2 Van der Waals 
complex to perform quantal close-coupling calculations and compute the index of refraction as a 
function of the projectile velocity. We obtain good agreement with the available experimental 
data. We show that the refractive index contains glory oscillations, but that they are damped by 
the averaging over the thermal motion of the N2 molecules. These oscillations appear at lower 
temperatures and projectile velocity. We also investigate the behavior of the refractive index at low 
temperature and low projectile velocity to show its dependence on the rotational state of N2, and 
discuss the advantage of using diatomic molecules as projectiles. 
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I. INTRODUCTION 

The refractive index is usually used to describe the 
propagation of light through a medium, but it can also 
be applied to the propagation of matter waves through a 
dilute gas. In that case, the refractive index is a complex 
quantity that depends on collisions between the parti- 
cles of the propagating wave and the gas The re- 
fractive index can be measured using atom or molecule 
interferometry and can help to determine properties of 
atom-atom or atom-molecule interaction potentials [li- 
lt is proportional to the complex forward scattering am- 
plitude f(k,0). The imaginary part describes the atten- 
uation of the beam and is related to the total scattering 
cross section, while the real part is associated with the 
phase shift due to the interaction between projectile par- 
ticles and atoms or molecules of the buffer gas. 

The first experiment on matter wave interferometry 
was realized with sodium atoms propagating in gas cells 
containing various gases Q . Subsequent experiments led 
to the observation of glory oscillations in the index of 
refraction as a function of the beam velocity for sodium 
matter waves in a medium of noble gases or nitrogen Q , 
and to the investigation of Na2 matter waves 0] . These 
works motivated studies on the theory of the refractive 
index for atoms which were used in combination 

with accurate interaction potentials to calculate the re- 
fractive index of various noble gases as a function of the 
sodium beam velocity [Iol - fl3j . A comparison between the 
experimental data and the theoretical calculations can be 
found in Ref. More recently, the refractive index of 
noble gases for lithium matter waves was measured (l5j 
and investigated theoretically [H, [l7| ■ 

Theoretical analysis of the matter wave refraction has 
been carried out mostly for atomic projectiles and bath 
gases. Collisions of molecular species may involve transi- 
tions between internal ro- vibrational degrees of freedom 



of a projectile or/and bath gas particles and the compu- 
tation of the refractive index becomes a significant prob- 
lem. In this work, we compute the index of refraction 
for sodium matter waves propagating through a gas of 
molecular nitrogen Nq. This system, for which experi- 
mental data are available over a wide range of projectile 
velocities has not been investigated theoretically. A 
comparison with the experiment also provides a check of 
the quality of the ab initio potential for the NaN2 com- 
plex. In Sec. mi we recall the definition of the refractive 
index in terms of the scattering amplitude, and we de- 
tail the calculation of the amplitude for atom-molecule 
collisions. In Sec. we compare our calculations to 
the experimental data and discuss the dependence of the 
scattering amplitude on the rotational level at low tem- 
perature and projectile velocity. 



II. THEORETICAL METHODS 

A. The refractive index 

In an interferometry experiment, a beam of projectile 
atoms with mass m p and velocity v p propagates through 
a cell containing a gas of target particles of density rit 
and mass m*. This gas can be described as an effective 
medium with a complex refractive index that depends 
on the scattering amplitude. Two expressions for the 
refractive index have been used in the literature. The 
first was derived by Champenois et al. Q, 
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while Forrey et al. wrote the refractive index as [11| 



n = l + 2Trn t —< — 



(1) 



(2) 



Present address: Service de Chimie Quantique et Photophysique, 
Universite Libre de Bruxelles, Belgium. 



In Eqs. © and ©, f(k r ) = f(k r ,9 = 0) is the forward 
scattering amplitude calculated in the center of mass 
frame, k r is the wave number for the relative motion 
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of the colliding particles, hk r = pv ri p is the reduced 
mass of the system, and k p is the wave number of the 
projectile particle in the laboratory frame, Kk p = m p v p . 
The brackets (•) denote a thermal averaging over the dis- 
tribution of target atom velocities. The real part of the 
refractive index is related to the phase shift accumulated 
by the matter wave due to scattering with the gas, while 
the imaginary part corresponds to the attenuation of the 
beam. All the results presented in this work were calcu- 
lated using the expression (0). Eqs. (TT|) and © can yield 
very different results at low projectile velocities, as was 
discussed by Champenois et al. [9J. However, for the sys- 
tem under consideration we found that the discrepancy 
between the two formulas was small even at low projectile 
velocities. 

While the real and imaginary parts of the refractive 
index can in some cases be measured accurately [15j . it 
is useful to define the quantity p as the ratio of the phase 
shift to the attenuation, 
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Re(n - 1) 
Im(n) 



(3) 



The main advantage of measuring p instead of n is that 
it is independent of the gas density Q • 



B. Thermal averaging 

While the projectile atoms have a definite velocity, this 
is not the case for the target atoms and a thermal aver- 
aging must be carried out. The thermal distribution for 
the velocity of the target atoms or molecules in the lab- 
oratory frame v t is given by a Maxwellian function: 



exp (-£ 2 v t 2 
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where T is the temperature of the target gas. In the 
experiment Q, T = 300 K. However, the scattering am- 
plitude is calculated in the center of mass frame, so that 
it is necessary to express this distribution as a function 
of the relative velocity v r = v p — v t . After integration 
over the angular part, the distribution (j4]) becomes 



P(v r , v p ) = -j= — exp {-p 2 {v 2 + v 2 )) sinh (2/3 2 v p v r ) 

V 77 v p 

(5) 

and the averaging of the scattering amplitude is given by 

(/(M> = / f{k r )P(v r ,v p )dv r (6) 
Jo 

C. Calculation of the scattering amplitude 

We calculated the forward scattering amplitude using 
the quantum close-coupling method [18| while treating 
N2 as a rigid rotor. The experimental temperature is 



much lower than the vibrational spacing of N2 so that 
all collisions occur in the ground state of vibrational mo- 
tion v = 0. The probability of excitation of N2 to higher 
vibrational levels is small at the experimental collision 
velocities and these channels may be ignored in calcula- 
tions of the refractive index. The rotational energies of 
N2 were obtained using the rotational constant [ljl [2(| 
B e = 1.99824 cm- 1 . 

In the close-coupling method, the wave function is ex- 
panded in terms of radial and rotational functions. This 
expansion is then inserted into the Schrodinger equation, 
leading to a set of second order equations in which the 
couplings are represented by matrix elements of the inter- 
molecular potential V(R, 9) [l8|, [2l[ . The close-coupled 
(CC) equations take a simple form if the total angular 
momentum J = j + 1 (where j is the angular momentum 
of N2 and 1 is the orbital momentum of the collision) is 
introduced. The total angular momentum is conserved 
during the collision, so that the coupled equations are 
block-diagonal in J. The S matrix elements S^,^ are 
obtained from the solution of the CC equations, and the 
scattering amplitude for the transition from an initial ro- 
tational state j with projection rrij to a final rotational 
state j' with projection my is given in terms of the S 
matrix elements as [l8| 

j I J 
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In this equation, Yv m , (9) are the spherical harmonics, M 
is the projection of J, kj = \j2p,E r jh = ^j2[i{E — tj)/h 
is the wavenumber in the entrance channel with energy 
ej, E r is the collision energy, and E is the total (kinetic 
plus rotor) energy. 

For the special case 9 = 0, 5W,,(0) = y^^,,o 
and the forward scattering amplitude becomes 

f jmj ^ fmj (E r , 0)-^EE i'"'V(2* + 1)(2/' + 1) 

x (-)'+''+'"+'" (2 J + l) ( 

In matter wave experiments, the detected particles 
have the same translational momentum as the projec- 
tiles. Therefore, we only need to compute the elastic 
scattering amplitude 



j I J 
rrij —M 
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(E r ,0), (9) 
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which is related to the total collision cross section via the 
optical theorem 



An 

Ptot,j(E r ) = — lm[f e i tj (E r ,Q)] 



(10) 



The dependence of the scattering amplitude on the ro- 
tational quantum number j must also be taken into ac- 
count. If we assume a thermal distribution of rotational 
levels, the amplitude entering Eq. ([2]) is given by 



f(k r ) = ^pj(T)/ clJ (fc r ) 



(11) 



where Pj(T) is the population of the rotational level j 
with energy ej 



(2j + l)exp(-6j/fcT) 
£,(2j' + l)ex P (-e//feT) 



(12) 



To obtain p(v p ) for a wide range of projectile velocities 
and temperatures, we computed the scattering amplitude 
for relative energies between 10 -8 and 10 4 cm -1 . At low 
energy, the CC equations were solved numerically [22} us- 
ing a log-derivative algorithm [23[ on a radial grid from 
^min = 3.5ao to i? max = 150ao for all total angular mo- 
menta J until convergence of the sum in Eq. ((5|). At 
high energy, we used the coupled states (CS) approxima- 
tion [24[ as implemented in the non-reactive scattering 
code MOLSCAT J2j|. The coupled states method leads 
to a drastic reduction in computational time and is par- 
ticularly useful at high energy. Above 1500 cm -1 , the 
agreement between the elastic cross sections calculated 
via the CC or CS methods was within 1%. 



D. Potential energy surface of NaN2 

We adopted the 2D potential energy surface (PES) cal- 
culated ab initio by Loreau et. al. and presented in 
Ref. (26|. The potential V(R,0) was obtained using the 
coupled-cluster method and a large basis set, with the N- 
N bond fixed to its equilibrium distance r e = 2.0743ar> 
The potential has a global minimum of 26.92 cm -1 at 
a distance R e = 10.47ao between Na and the center of 
mass of N2 and at an angle 9 e = 45.0°. The poten- 
tial has a saddle point for 9 = 90° which lies about 4 
cm -1 above the global minimum. The long range part 
of the potential was constructed assuming a Cq/ R e form 
and using calculated dispersion coefficients [26j , and was 
joined smoothly to the ab initio points using a switch- 
ing function. To perform scattering calculations, this 
potential was expanded in a basis of Legendre polyno- 
mials, V(R, 9) = J2\ V\{R)P\ (cos 9), retaining terms up 
to A = 12. Additional terms did not affect the cross 
sections or scattering amplitudes. 



III. RESULTS AND DISCUSSION 
A. Comparison with the experiment 

The parameter p(v p ), defined by Eq. ©, is shown 
in Fig. [T] for projectile velocities between 800 and 3000 
ms _1 , together with the experimental values from Ref. 

. These experimental data were obtained using 200 nm 
gratings at low velocities and 100 nm gratings at higher 
velocities. We obtain very good agreement between our 
calculations and experiment for projectile velocities up to 
around 2000 ms _1 . However, the theoretical values of p 
are lower for velocities higher than 2500 ins -1 . We have 
no explanation for the discrepancy but our calculations 
show that it is not due to the contribution of inelastic 
rotational transitions. It is also unlikely to be due to 
uncertainties in the calculated potential energy surface, 
despite the sensitivity of p to it [llj . Discrepancies also 
occur in the simpler case of sodium matter waves travel- 
ing through noble gases [14] . 
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FIG. 1. Comparison between our calculations (full line) and 
the experimental values (using 200 nm or 100 nm gratings [J]) 
of p as a function of the projectile velocity v p . The tempera- 
ture is T = 300K. 

In contrast to the noble gases, the experimental data 
and our calculations for N2 show no sign of glory oscilla- 
tions. The absence of glory oscillations in the refractive 
index despite their presence in the Na-N2 elastic collision 
cross section [26| occurs because they are suppressed by 
the averaging over the thermal distribution, as found for 
other systems [l(| El, E3] • The averaging depends on the 
temperature as well as on the mass of the target, and 
the damping would be weaker for heavier particles or for 
lower temperatures. However, even with a temperature 
lower than 300 K, glory oscillations still would not be 
observable for the range of velocities used in the experi- 
ment. This is illustrated in Fig. [21 which displays p(v p ) 
for several temperatures and for v p between and 3000 
ms _1 . It can be seen that p does not present any glory os- 
cillations for projectile velocities between 1000 and 3000 
ms _1 , even at low temperature. However, for v p below 
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1000 ms , the oscillations start to appear and become 
more marked as the temperature decreases. 
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FIG. 2. Behavior of the ratio p(v p ) for various temperatures. 



This can be understood by looking directly at the for- 
ward scattering amplitude, shown in the top panel of 
Fig. |3] for energies between 0.5 and 10 4 cm' 1 . Glory 
oscillations are indeed present in both the real and imag- 
inary parts of the amplitude, but they are restricted to 
low energy. The scattering amplitude is shown for the 
rotational level j = 7, which is the most populated at 
T = 300 K. The amplitude is almost independent of the 
rotational level over the energy range considered, as al- 
ready noted for the scattering cross sections [26(. The 
j— averaging has therefore no effect on the value of 
p under the experimental conditions. 

The distribution ([5]) is shown in the bottom panel of 
Fig. |3]for the lowest and highest projectile velocities used 
in the experiment. These distributions probe the scatter- 
ing amplitude in the region where no glory oscillations are 
present, so that the averaging of f(k r ) does not yield any 
oscillations in p(v p ). If we decrease T and v p , it becomes 
possible to scan the glory oscillations that are present at 
lower collision energy. At the lowest temperature shown 
in Fig. [U (T = 2K), the thermal distribution flSJl is so 
narrow that the features of the scattering amplitude are 
directly reflected in p(v p ); in particular, each maximum 
(minimum) of p corresponds to a maximum (minimum) 
in the scattering amplitude. 

The scattering amplitude also shows a resonance struc- 
ture at energies below the depth of the potential (27 
cm -1 ), which has been discussed in Ref. [2g]. This struc- 
ture depends on the rotational level j. 

Finally, we also calculated p(v p ) using Eq. ([T]) instead 
of Eq. © . We found that the two formulas give results 
that agree within a few percent, even at low projectile 
velocities. 
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FIG. 3. Top panel: real and imaginary parts of the forward 
scattering amplitude f(k) as a function of the relative en- 
ergy for j = 7. Bottom panel: distribution function (5| as a 
function of the relative energy for T = 300 K. 



B. Molecular matter waves 

An interesting question for matter wave interferome- 
try involving molecules is whether exchanging the pro- 
jectile and target (i.e., using a beam of molecules trav- 
eling through a cell of atoms) can provide new insights, 
in particular related to the determination of properties 
of the interaction potential. For example, it has already 
been pointed out in the case of atom interferometry that 
the magnitude and sign of the scattering length can be 
obtained by measuring p 0, [l2j]. Apart from the dif- 
ference in mass, which modifies the thermal distribution 
([5]), the most obvious consequence of using molecules as 
projectiles is that if the molecules can be prepared in 
a definite rovibrational level, then an interferometry ex- 
periment could lead to insights into the nature of the 
interaction potential involving those specific levels. This 
would be particularly interesting at low temperature as 
p is then expected to be strongly dependent on v and j. 

In the following, we will use the scattering amplitudes 
calculated in the previous section to investigate the be- 
havior of p at low projectile velocity and temperature. 
We assume that a beam of nitrogen molecules prepared 
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in a given rotational state travels through a gas cell of 
sodium atoms at a temperature T = 10~ 8 K, so as to 
avoid any effect of the thermal averaging. As a conse- 
quence, the two formulas ([TJ and @ give the same value 
for p. 

In the s-wave regime, the scattering amplitude can be 
written in term of the phase shift So as 



f(kr) 



exp(2iS ) — 1 
2z ky 



So 

k, 



(13) 



since So = —k r a in the limit k r — ¥ 0, where a is the 
scattering length. In the case of atom-molecule collisions, 
the molecule may be in a ro-vibrationally excited state 
and the scattering length can be complex, a — a ~ i/3. 
The imaginary part j3 is related to the inelastic cross 
section by f3 = fc r 0mci/47r in the limit k r — > [27j . and 
the elastic cross section is given by <r c \ — An (a 2 + f3 2 ). 
Using Eq. (fT3j) , we see that the s-wave amplitude is given 
in the limit k r — > by / = —a + i/3, so that p — —a//3. 

In the case of Na-N2 scattering, (3 = for j = and 
j = 1 so for these levels p will tend to ±oo when k r — > 0. 
On the other hand, for j > 1 quenching to levels j' < j 
is possible and p will be constant. Its value will depend 
on the relative contribution of elastic and inelastic scat- 
tering to the total cross section. We computed the ratio 
p at low velocity and temperature for various rotational 
levels. The results for j = 0, 3, 5 and 6 are illustrated in 
Fig. SJ For j — 0, a is positive and p decreases to — oo, 
as expected. For j = 3 or j = 5, a is negative and has 
a similar magnitude as f3 so that p w 1, while for j = 6 
a is positive and much greater than f3. The behavior of 
p can be understood from Eq. (1131) by retaining higher 
order terms in the expansion of exp(2i<5o). For exam- 
ple, keeping second order terms in the fc-expansion of the 
scattering amplitude, we find that the real part of the 
scattering amplitude for j = is still Re(/) = —a, while 
the imaginary part is Im(/) = k r a 2 . Therefore, for this 
level p can be approximated at low energy by [281 ] 



1 

k r a 



(14) 



This function is shown in dashed line in Fig. |4]and agrees 
very well with the full calculation up to v p ~ 1 ms -1 . 

It should also be noted that at projectile velocities 
v p > 20 ms _1 , the value of p becomes almost independent 
of the rotational level j, as indicated previously. There- 
fore, it is clear that at high temperatures and projectile 
velocities, the Na-N2 system does not represent a good 



experimental case to extract information about the inter- 
action potential from matter wave interferometry. Other 
systems such as Na2-He (recently studied by Bjorgen et 
al. [13]) could be be more promising as p seems to have 
a dependence on the rotational level at higher projectile 
velocities. 
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FIG. 4. Behavior of p(v p ) for various rotational levels at 
low projectile velocity and a temperature T = 10 -8 K. The 
dashed line corresponds to Eq. (TT4]) . 



IV. CONCLUSION 

We computed the index of refraction of sodium mat- 
ter waves propagating through a molecular nitrogen gas. 
We obtained the scattering amplitude by performing fully 
quantal calculations using an accurate ab initio potential 
energy surface for the NaN2 complex. The value of the 
ratio p of the real to the imaginary part of the scatter- 
ing amplitude agrees well with the experiment except at 
the highest experimental projectile velocities. We showed 
that no glory oscillations were observed in the experiment 
due to the fact that both the temperature and the sodium 
velocity were too high. Finally, we discussed the behav- 
ior of p in the ultracold regime and its dependence on 
the rotational level of N2 . We showed that for ultracold 
collisions the refractive index is very sensitive to the val- 
ues of the molecular rotational quantum numbers. This 
sensitivity may be used for an accurate determination of 
parameters of the j— dependent interaction potentials in 
atom-molecule systems. 
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